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For the one-dimensional Holstein model, we show that the relations among the scaling exponents 
of various correlation functions of the Tomonaga Luttinger liquid (LL), while valid in the thermo- 
dynamic limit, are significantly modified by finite size corrections. We obtain analytical expressions 
for these corrections and find that they decrease very slowly with increasing system size. The inter- 
pretation of numerical data on finite size lattices in terms of LL theory must therefore take these 
corrections into account. As an important example, we re-examine the proposed metallic phase 
of the zero-temperature, half-filled one-dimensional Holstein model without employing the LL rela- 
tions. In particular, using quantum Monte Carlo calculations, we study the competition between 
the singlet pairing and charge ordering. Our results do not support the existence of a dominant 
singlet pairing state. 

PACS numbers: 71.10.Fd, 71.30.+h, 71.45.Lr 
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I. INTRODUCTION 

The study of quantum fluctuations arising from dy- 
namical phonons and their effects on the charge order- 
ing instability of one-dimensional systems has a long 
history^ One of the simplest models to describe a one- 
dimensional solid under the influence of the quantum 
fluctuations is the molecular-crystal or Holstein model^i 3 - 
which is given by the following Hamiltonian: 



where 



H — H e i + H p h + H e i_ p h, 



H e i = -t^2(c\ s Ci + i >s + H.c.) , 

i,s 



Hei-ph — .9cp y 



(1) 

(2) 

(3) 
(4) 



In this model, the electrons move in a tight-binding 
band interacting with the lattice via the coupling between 
the local electron density and the local lattice displace- 
ment. In the Hamiltonian, c\ s (cj iS ) creates (annihilates) 
an electron at site i with spin s, and n^ s is the electron 
number operator. The ions are described by momentum 
Pi and displacement for site i. The mass of the ions 
is M, and K is the harmonic lattice coupling constant, 
so that the phonon frequency, ujq — \J K/M. Following 
convention, we shall call these dispersionless optical ex- 
citations "Einstein phonons" . The hopping energy of the 
electrons between nearest neighbors is t, and g cp is the 
electron-phonon coupling energy. In what follows we will 
use units such that t = M = 1, except where otherwise 
specified. For the half-filled band, early work suggested 
that the Holstein model displays dimerization and that 



long range charge ordering sets in at any finite phonon 
frequency.— ~— These conclusions have been challenged by 
various studies which find that for small enough electron- 
phonon coupling, the system has zero charge gap and 
thus no long range order.— ~— 

A recent quantum Monte Carlo study2^ further sug- 
gested that the dominant power-law ordering is the sin- 
glet BCS superconducting pairing (SS) correlation. This 
conclusion was based on data from finite size simulations 
which shows that the Luttinger charge exponent (K p ) is 
larger than one^ (It is important to realize at the out- 
set that this is not an off-diagonal long range ordering 2 -^, 
as breaking the U(l) gauge invarianc o 28 ' 29 to generate a 
mass gap is forbidden at one dimension.— ~— SS should 
thus be understood as a charge gapless, dominant off- 
diagonal power-law decaying ordering.) In the (possi- 
bly singular) limit of infinite phonon frequency, the Hol- 
stein model becomes equivalent to the attractive Hub- 
bard models For the attractive Hubbard model at half- 
filling, the SS and the charge density wave (CDW) or- 
der parameters are degenerate because of the particle- 
hole symmetry. Moreover, both order parameters possess 
continuous symmetry with respect to the spin rotation. 4 
Therefore, long range ordering of either sort is forbidden 
in this case. However, from the numerical simulations for 
finite size lattices, K p is always larger than one and con- 
verges to one very slowly with increasing system size ^2 
The similar behaviors of K p for both the attractive Hub- 
bard model and the Holstein model at finite phonon fre- 
quency have led to the suggestion that the Holstein model 
even at finite phonon frequency may have a charge gap- 
less metallic phase with degenerate SS and CDW correla- 
tions, even though K p > 1 from finite size calculations. 33 
We note that this coexistence has been proposed years 
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In the present article we present analytical and numer- 
ical calculations of the LL relations in the ID Holstein 
model with finite size corrections. Our results extend 
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to Einstein phonons previous studies involving acoustic 
phonon o 35 ' 36 and complement prior multiscale functional 
renormalization group studies of retardation effects due 
to Einstein phonons M- In essence, we seek to answer three 
important questions raised by previous studies: (I) Are 
the LL scaling relations valid in the ID Holstein model? 

(II) What are the dominant correlations in this model? 

(III) What can be said about the possibility of a charge 
gapless metallic phase at half-filling for weak electron- 
phonon coupling? As we shall see, we are able to provide 
clear answers to the first two questions but can only com- 
ment on the third. 

In section II of the paper, we begin our study by de- 
scribing the forward scatterings of the Holstein model in 
terms of bosonic fields. In section III, we solve analyt- 
ically the forward scattering part of the Holstein model 
to show that electron-phonon coupling modifies the rela- 
tions between scaling exponents , and we calculate these 
modifications for finite size systems. This enables us to 
answer question (I). In section IV, we employ quantum 
Monte Carlo calculations for the full Hamiltonian (Eq. 
[I} to resolve Question (II) and to gain some insight into 
Question (III). In section V, we discuss the reasons of 
why all the perturbative renormalization group studies 
obtain long range CDW phase and under what circum- 
stances the metallic phase can be obtained. Section VI 
contains a brief conclusion. 



II. FORWARD SCATTERINGS IN THE ID 
HOLSTEIN MODEL 

To answer the question about the validity of the LL 
relations for the Holstein model in the gapless regime, 
including the gapless half-filled case if present, it is suffi- 
cient to diagonalize exactly the forward scatterings ( scat- 
terings with zero momentum transfer, often denoted as 
<72 and <?4 in the "g-ology" terminology)^ of the Hamil- 
tonian. The neglect of anything beyond forward scat- 
terings is unquestionably an assumption which cannot 
be a priori justified in the full model, but this neglect 
is the fundamental conjecture on the description of the 
gapless one-dimensional system: any contributions be- 
yond forward scatterings can be incorporated into the 
renormalization of the forward scatterings. Therefore, 
although we are considering forward scatterings only, it 
is completely general for the description of any gapless 
regime regardless of the filling. In what follows, the right 
hand side of the symbol, sa, is understood to include only 
those terms in which the forward scatterings exclusively 
are retained. 

By formulating the one-dimensional problem in the 
continuum limit and retaining only the low-energy de- 
grees of freedom, H e i , can be expressed in term of bosonic 
fields, tp u and their corresponding momenta n„ with ex- 
plicit spin(cr) - charge(p) separation! 39 ' 40 



(5) 



= Jdx[r 



For the Holstein model at half filling, the Fermi points 



are at kp = —tt/2 and tt/2. Therefore u a 



2t, and 



K a = K p = li 35 i 36 ' 39 The values of these parameters, 
u a ,u p , K a , K p , can be renormalized by the interactions. 
However, as long as the LL description is valid, the scal- 
ings of the various correlation functions are determined 
solely by these parameters. 39 ' 40 In particular, the CDW 
correlation and the SS correlation scale as Ocdw{x) oc 

X -K p -K a _ x «cdw an( J SS {x) (X X~ K ^~ K ° = X ass , 

respectively. This implies acDW&SS = 1 ; when a spin 
gap appears, in which case K a = 0. 

Consider the electron-phonon coupling. As we include 
only forward scatterings, only the phonons with momen- 
tum k = are coupled with the electrons. The Hamilto- 
nian for the phonons with momentum k = is 



H, 



1 



ph 



dx[c _1 n 2 



''0 



(6) 



where £ is the mass density, do is the phonon field with 
momentum k = 0, and is the conjugate momentum. 
The Hamiltonian for the electron-phonon coupling is 



Hei-ph ~ 9c P ^2 J dx(i[>lip s d ). 



(7) 



,35,36.39 



By keeping only the term which is linear in the 
bosonic field to express the density operator, ipli/j s , we 
obta 



(8) 



Hei-ph ~ g cp / dx\ -{d x (pp)d 



The remaining terms will generate the back scatterings 
and Umklapp scatterings. They also lead to spin-charge 
coupling due to the retardation^ 

When only forward scatterings are considered, the 
electron-phonon coupling affects only the charge part of 
the Hamiltonian, and we will thus focus on the charge 
part only in the following. The spin part is unchanged 
from the standard Luttinger mode L 35 ' 36 ' 39 The charge 
part of the Hamiltonian (Eq. [T]) arising solely from for- 
ward scatterings is 

H forward ~ J dx[u p K p Il 2 p + ^-(d x ip p ) 2 } + (9) 

i J dxlC 1 !! 2 ^ + C^o d o] +5o P J dx^(d x ip p )d . 

Since the phonon fields appear at quadratic order, we 
follow the work by Martin and Loss ; 35 ' 36 who considered 
the case of acoustic phonons (as opposed to the Einstein 
phonon case we are considering here) to integrate out the 
phonon fields, d$, which results in the collective propa- 
gator for the charge field, ip p , 



D p (k,u) 



1 



KpUp 



u 2 p k 2 ) 



b(koj u p ) 



(10) 
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where b = 2K p g^ p /nu p ujQ. The first part of the prop- 
agator is the contribution from the bosonized electrons, 
and the second part comes from the electron-phonon cou- 
pling. The crucial difference between Einstein phonons 
and acoustic phonons is that the coupling between the 
two bosonic modes (tp p and do) is symmetric for acous- 
tic phonons (d x ip pd x d{f)r^^, but asymmetric for Einstein 
phonons (d x ip p do). As a result, the dispersion relations 
for the collective modes are linear for acoustic phonons 
but nonlinear for Einstein phonons. Therefore, the acous- 
tic phonon case can be solved exactly, but the Einstein 
phonon case, in general, cannot be solved exactly. 



III. MODIFICATIONS OF SCALING 
EXPONENTS 

It is generally believed that the spin part of the half- 
filled Holstein model is always gapped due to the back- 
scatterings ^ 25 ' 33 (The spin-gapped but charge-gapless 
LL is often referred as Luther-Emery liquid4^) As a con- 
sequence, the Luttinger spin exponent K a renormalizes 
to zero. We take this as an assumption and will not elab- 
orate on it here. Therefore, the boundary of the transi- 
tion between the charge ordering and the singlet pairing 
is determined by the points where the Luttinger charge 
exponent crosses the value one, providing that the Lut- 
tinger liquid description survives. 

It has previously been shown that the LL scaling 
relations are modified by the coupling with acoustic 
phononsi 35 ' 36 The analogous exact relations cannot be 
obtained for the Holstein model, which has Einstein 
phonons. Nevertheless, we can still extract the power-law 
decaying terms asymptotically (i.e., long-distance, low- 
energy limit, x — > oo, k — > 0). We find 



&CDW 



ass 



-K, 



(11) 
(12) 



to the lowest order in the momentum. At this order, the 
velocity of the lower branch, u_ is u p \J\ — b. Hence the 
relation acDwass = 1 is preserved^ We reiterate that 
these asymptotic results are strictly correct only at long- 
distance, even apart from any oscillatory contributions. 
In addition, the power-law decay of the correlations is 
also valid at long distances only. 

The exact results for the charge density, Ocdw, and 
the singlet pairing, Ossi correlation functions can be ex- 
pressed in term of two series: 



-1 



\n{0 C Dw{b,x)) = — 



2tt 



k 



1 



HOss{b,x)) 



[1 — cos(kx) exp(- 

-y 

2tt ^ 

k 



K p uj- 



■w_t)], 
(1-6V 2 



(13) 

■ ,2 



where 



<4(*)= 



luI + {upkf ± y/[u* - {upk) 2 ] 2 + Ab(io u p ky 



(15) 



We have omitted oscillatory contributions. In addition, 
the upper branch of the collective modes is ignored due to 
the gap with magnitude, this differs from the acoustic 
in which both branches are gapless and must be 



casi 



,35.36 



included. 

To investigate the validity of the relation acDW a ss — 
1 at finite distance we take the continuum limit, which 
amounts to replacing the summation over the momentum 
by an integral; and we have to introduce the low momen- 
tum cutoff c 35 ' 36 and the high momentum cutoff is set 
to 2u p . For the following numerical calculations, we fix 
the phonon frequency, ljq = 1. The asymptotic results 
indicate that the relation, acDwass = 1, is valid inde- 
pendent of any parameter. In particular, it is valid in- 
dependent of the electron-phonon coupling, which enters 
via the parameter b. Therefore, the ratio of the product 
of the singlet pairing and the charge ordering correlations 
at b ^ and b = 0, 



R{b,x) 



ln[Q C Dw{b, x)]ln[0 S s(b, x)} 
ln[O C Dw{0,x)]ln[O S s(0,x)]' 



should be equal to one when x — » oo, that is, 

Um R{h x) = WMQ^jMO) = L 

z->oo acDw{b — 0)ass{b = 0) 



(16) 



(17) 



cos(kx) exp(— ui—t)], 



(14) 



While Eq. [17] is exact in the asymptotic limit, using 
the analytical expressions we obtained for Oss{b,x) and 
OcDw(b,x), we find that R{b,x) deviates from one for 
non-zero scaled electron-phonon couplings, b; even at 
rather long distances, x. This is shown in detail in Fig. 
Q] Importantly, this deviation depends not only on the 
distance but also on the value of the electron-phonon cou- 
pling. The ratio is clearly larger than one and seems to 
converge to the asymptotic value, one, very slowly and 
at distances that are likely beyond the reach of exist- 
ing numerical calculations, which typically study systems 
up to a few hundred sites. This behavior is also likely 
to preclude a reliable scaling estimate of its asymptotic 
value. Our results suggest that the product of the expo- 
nents of the singlet pairing and of the charge ordering is 
larger than one at finite distance for non-zero electron- 
phonon coupling. Consequently, observing a charge ex- 
ponent larger than one does not automatically imply that 
the singlet pairing correlation is the dominant instability. 
Thus although acDwass = 1 in the asymptotic limit, 
in a finite size system the value of acDwass n ot only 
depends on the system size but also on the forward scat- 
terings strengths. 

The deviation from one is enhanced by the scaled 
electron-phonon coupling, b, and suppressed by the dis- 
tance, x. These results are consistent with the charge 
exponent obtained in a Quantum Monte Carlo study, 25 
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FTO 1 • Pint of the ratio R(h r\ = ln lOcDw(b,x)]lnlO ss (b,x)] 
versus the scaled electron-phonon coupling, b, at various dis- 
tances, x, for t — > + 



in which the charge exponent is larger than one and in- 
creases with electron-phonon coupling. In addition, the 
Quantum Monte Carlo calculations also found that the 
charge exponents decrease with increasing system size^ 

Our results show that even if all the non-linear op- 
erators, which are always present and especially impor- 
tant in the half-filling case, are irrelevant and renor- 
malized to zero and the Holstein model is in the gap- 
less regime which can be described in terms of LL, 
the charge exponent alone at finite distance cannot re- 
solve the competition between the charge ordering and 
the singlet superconducting pairing. In fact, the "two- 
cutoff' renormalization group^2r— calculations have 
shown that the charge gap is always opened by the 
Umklapp scatterings^ Similarly, multiscale renormal- 
ization group approaches which systematically renormal- 
ize the frequency dependences of the couplings^— also 
ruled out the existence of dominant singlet pairing due 
to the contribution from the dynamical Umklapp scat- 
terings at high energy^i Further, a recent study once 
again confirmed the CDW phase and further discussed 
the quantum-classical crossover, 50 

From the results of the correlation functions for the 
forward scatterings of Eq. [TJ it is quite clear that using 
the relation, otcDW&ss = 1, to decide the dominating 
correlations in the gapless phase of the half-filled Hol- 
stein model can lead to misleading results. We reiterate 
that the violation of this relation at finite system sizes is 
not because of the oscillating terms from normal order- 
ing but because of the non-linear collective modes from 
the asymmetric coupling between the electrons and the 
dynamical phonons. 



IV. MONTE CARLO CALCULATIONS OF THE 
CORRELATION FUNCTIONS 

Our previous results make clear that we cannot rely 
on the LL relations for the thermodynamic limit to de- 
termine the correlation exponents from finite size sim- 
ulations of the Holstein model. Hence we turn instead 
to a numerical study of the full model (Eq. 1) and use 
the projector determinantal Monte Carlo 5 -^— to deter- 
mine the "ground state" and its correlation functions of 
finite size systems. The "ground state" (|GS)) is ob- 
tained by projecting the free fermion trial state (I^t)) 
as |GS) = lim ( 3_ ! . 00 exp(— (3H)\^^), in practice (3 is set 
between 10 — 20. We use the |GS) to calculate the pair- 
ing correlation, Oss, and the charge ordering correlation, 
Ocdw, in the real space directly, 



N 



CDW (x) = hj2( ni n i+x -l)\; (18) 

i=l 
1 N 



N' 

i=l 

where N is the number of lattice sites. In the Monte 
Carlo calculation, we set t = 1, M = 1/2, and ojq = 1, 
so that the g ep defined in Eq. [T] is the same as the 
electron-phonon coupling g defined in the Ref. [2o| . 
We show the correlation functions at different electron- 
phonon couplings in Fig. [2] At zero electron-phonon 
coupling, the Holstein model is a free fermion theory, 
and the singlet pairing and charge ordering are degen- 
erate (the 1st row in the Fig. [2]). When the electron- 
phonon coupling is turned on, the charge ordering tends 
to decay more slowly, while the singlet pairing tends to 
decay faster, than in the free fermion case. The differ- 
ence between them is expectedly small, possibly due to 
the exponentially small charge gap with respect to the 
electron-phonon coupling A Turning to larger electron- 
phonon couplings, the difference between the two differ- 
ent orderings becomes more obvious. As we are using 
systems with modest sizes (22 and 46 sites with peri- 
odic boundary condition) only, it is difficult to determine 
whether the charge ordering is long range or decays as a 
power law, but the data clearly show that the charge or- 
dering is enhanced and the singlet pairing is suppressed 
as the electron-phonon coupling is increased. For fairly 
large electron-phonon couplings, e.g., g cp = 1.0 (the 6th 
row in Fig. [5]), our data show a clearer signal that the 
charge ordering tends to become long range, as the charge 
correlation changes only very slowly beyond a few lattice 
sites. 

We also calculate the Luttinger charge exponent by the 
same method used previously for the Holstein model, 25 
which found K p > 1. Using the set of data obtained in 
the same run for the correlation functions presented in 
Fig. [21 we reproduce the result that the Luttinger charge 
exponent is clearly above one at least for g cp < 0.6 (see 
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FIG. 2: Singlet pairing correlations (open circles) and charge 
ordering correlations (open squares) functions as a function 
of real space distance for g cp = 0.0 (the 1st row), g ep = 0.2 
(the 2nd row), g cp — 0.4 (the 3rd row), g cp = 0.6 (the 4th 
row), g cp = 0.8 (the 5th row) and g cp — 1.0 (the 6th row). 
The left and right columns are for the lattices of N = 22 sites 
and N = 46 sites, respectively. 



Fig. [3]). Using the LL relations would imply the dom- 
inant singlet pairing state for the model in this region; 
however, from the correlation functions we know that the 
singlet pairing decays faster than the charge ordering and 
is thus not dominant. The K p > 1 obtained from finite 
size calculations is not incompatible with the CDW cor- 
relation decaying slower than that of the SS. 




FIG. 3: The Luttinger charge exponents K p and spin expo- 
nents K a for different electron-phonon couplings. iV is the 
number of lattice sites. 

From our Monte Carlo data, we can conclude that sin- 
glet pairing is not the dominant instability for any range 
of electron-phonon coupling (although, strictly speaking, 
numerical data can hardly provide a conclusive answer 
for g ep — > + ). This provides a reasonably convincing an- 
swer to Question (II) . The most prominent open question 
is whether the long range CDW ordering or equivalently 
the charge gap can be destroyed by the quantum fluctua- 
tions from the dynamical phonons as the electron-phonon 
coupling is decreased below a certain value. 



V. DISCUSSION 

Perturbative renormalization group approaches have 
consistently shown that the charge part is always gapped 
and thus the CDW is long range ,l & 37 i 44 1 46 £ Thig ig ulti _ 

mately related to the fact that the continuous spin rota- 
tion symmetry is broken by any finite phonon frequency 
which leads to a relevant Umklapp scattering. This re- 
sult already appears in the perturbative renormalization 
group approaches within a second-order perturbation ap- 
proximation. For a metallic phase to be obtained in 
a perturbative renormalization group approach for the 
Holstein model, in particular for degenerate CDW and 
SS phase to be obtained, the summation of all the higher 
order terms in the perturbation series would have to can- 
cel exactly the second order contribution to restore the 
spin rotation symmetry. This is an unlikely scenario, and 
there is no convincing argument for why this should hap- 
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pen. 

On the other hand, to determine precisely the na- 
ture of the ground state from numerical calculations 
for very weak couplings is challenging, as the charge 
gap is exponentially small with respect to the electron- 
phonon coupling, even within the mean field approxima- 
tion (A p ~ exp(— l/<7g p )). Quantum fluctuations from 
the dynamical phonons will further suppress it, espe- 
cially at large phonon frequency^i^i^ Nonetheless, recent 
density matrix renormalization group studies did argue 
for the metallic phase, and further claimed that there 
are two types of metallic phases in both Holstein and 
Holstein- Hubbard modeL 21 i 22 One of the metallic phases 
is gapless in both the spin and charge sectors, which nor- 
mally implies that either triplet superconducting pairing 
or spin density wave ordering is the dominant power-law 
decaying ordering^, even though spin fluctuations are 
completely absent in the Holstein model. More elabo- 
rate calculations using Monte Carlo algorithms which can 
directly extract the pairing correlations for much larger 
system sizes may provide further insight.— 

Although we cannot rule out the possibility of a metal- 
lic phase from this study, it is quite clear that the Holstein 
model in one dimension does not support dominant su- 
perconducting correlations for the half-filled case. This is 
due to the relevant Umklapp scattering from the perspec- 
tive of the perturbative renormalization group t 6 ' 7 ' 37 i 50 
which is supposed to be reliable from weak to intermedi- 
ate coupling. In general, breaking the particle-hole sym- 
metry alone,— by, for example, including non-linear term 
for the phonon, is probably not sufficient to induce dom- 
inant superconducting correlation^ The crucial point is 
to weaken the contribution from the particle-hole chan- 
nel by removing the nesting condition. Theoretically, this 
can be achieved by introducing long range hopping terms 



or doping away from the half-filling 



VI. CONCLUSION 

In this paper we have examined the nature and ap- 
plicability of using the LL relations to determine the 
possible existence and the nature of a proposed metallic 
ground state of the one-dimensional half-filled Holstein 
model. We conclude by reiterating our answers to the 
three questions posed in our introduction. (I) The LL 
scaling relation (xcdwoiss = 1 is valid only in the limit 
of infinite size systems; for finite size systems there are 
corrections that lead to K p > 1. (II) For a wide range of 
electron-phonon coupling, the dominant correlation for 
finite size systems in the Holstein model is the CDW. 
(Ill) Although the CDW seems to be always the domi- 
nant correlation, without a direct study of the gap-which 
is expected to be exponentially small with respect to the 
electron-phonon coupling and thus very challenging to 
study numerically-we cannot rule out the possibility of a 
charge gapless "metallic" phase for weak electron-phonon 
coupling. 
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